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Abstract 

The mitogen-activation protein kinase ERK2 is tightly regulated by multiple phosphatases, including those of the kinase 
interaction motif (KIM) PTP family (STEP, PTPSL and HePTP). Here, we use small angle X-ray scattering (SAXS) and isothermal 
titration calorimetry (ITC) to show that the ERK2:STEP complex is compact and that residues outside the canonical KIM motif 
of STEP contribute to ERK2 binding. Furthermore, we analyzed the interaction of PTPSL with ERK2 showing that residues 
outside of the canonical KIM motif also contribute to ERK2 binding. The integration of this work with previous studies 
provides a quantitative and structural map of how the members of a single family of regulators, the KIM-PTPs, differentially 
interact with their corresponding MAPKs, ERK2 and p38oc. 
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Introduction 

The mitogen-activation protein kinases (MAPKs; ERK, p38 
and JNK) are cytosolic serine/threonine-specific kinases that 
coordinate the cellular response to a range of extracellular stimuli. 
Their activity is tightly regulated by the coordinated action of 
activating kinases (the MAP kinase kinases), deactivating phos- 
phatases (i.e., the kinase interaction motif protein tyrosine 
phosphatases or KIM-PTPs, the dual specificity MAPK phospha- 
tases or MKPs, among others) and scaffolding/ targeting proteins, 
which determine their subcellular location [1]. Once activated, the 
MAPKs phosphorylate a variety of substrates, which ultimately 
determine the cellular response [2]. Consequently, disruptions in 
MAPK signaling are correlated with multiple diseases, including 
cancer, autoimmunity and Alzheimer's disease, among others [3]. 

MAPKs have a bi-lobed architecture, with a five-stranded (3- 
sheet in the N-terminal lobe and six a-helices in the C-terminal 
lobe [4,5]. During the last two decades, cellular and structural 
studies have revealed that MAPKs contain a conserved docking 
groove, characterized by two distinct subsites (one charged, "Pchgj 
and the second hydrophobic, <I>hyd) [6-8]. Most regulators and 
substrates bind their cognate MAPKs at this groove via linear 
peptides, known as KIM- or D-site motifs (Kinase Interaction 
Motif), with the consensus sequence (R/K) 2 _ 3 -X 2 _ 6 -<I> A -X-<I> B (<t is 
any hydrophobic residue) [9-11]. One family of regulators that 
bind MAPKs via a KIM is the KIM-PTPs (Kinase Interaction 
Motif - Protein Tyrosine Phosphatases) [1 1-13]. Members of this 



family include STEP (PTPN5, striatal enriched protein tyrosine 
phosphatase), PTPSL/PTPBR7 (PTPRR, PTP STEP-like) and 
HePTP (PTPN7, hematopoietic protein tyrosine phosphatase). 
They regulate ERK2 (extracellular signal-regulated kinase 2) and 
p38oc by selectively dephosphorylating their activation loop 
tyrosines (active state complexes) and also by sequestering the 
kinases in the cytosol (resting state complexes). 

Over the last years different groups, including ours, have studied 
the interaction of MAPK regulatory proteins, as well as substrates, 
with different MAPKs, especially ERK2, p38 and JNK [12-17]. 
However, while a plethora of cell biological data is available, much 
less systematic interaction data using quantitative methods, such as 
isothermal titration calorimetry (ITC), has been reported. This is 
of special interest, as recent studies have shown that while KIM- 
peptides derived from KIM-containing regulatory proteins and 
substrates readily discriminate between ERK2/p38 and JNK, they 
discriminate poorly between ERK2 and p38 [7]. We recently 
reported detailed ITC, NMR as well as small angle X-ray 
scattering (SAXS) data of the interactions between the MAPK 
p38ot and the KIM-PTP family members [14,15]. Furthermore, 
we also compared the interaction the immune-specific KIM-PTP, 
HePTP, with p38c< and ERK2 [16,17]. These studies showed that 
HePTP (residues 15-339) as well as the HePTP KIM-peptide 
(residues 15-31) consistently bind more tightly to ERK2 (~7-times 
and ~4-times, respectively) than p38ot. To better understand the 
molecular determinants that allow regulatory proteins to discrim- 
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inate between ERK2 and p38oc, we investigated the interaction of 
PTPSL and STEP with ERK2 and then compared them with our 
previous study examining the same interactions with p38oc. Our 
data shows that both PTPSL and STEP also bind more strongly to 
ERK2. Furthermore, we find that within the family both STEP 
and PTPSL interact with ERK2 and p38a more strongly than 
does HePTP. To determine if this discrimination between p38 and 
ERK2 is specific for this regulatory protein family, we also 
compared the interaction of a MAP Kinase Phosphatase (MKP; 
DUSP16) MAPK binding domain with p38a [18] and ERK2. In 
contrast to what we observed for the KIM-PTPs, DUSP16 binds 
both MAPKs with similar affinities, demonstrating that other 
families of MAPK regulator proteins likely have distinct discrim- 
ination factors between the MAPKs. Finally, we correlated these 
differences with the low-resolution structures of the complexes 
using small angle x-ray scattering (SAXS) experiments. Collec- 
tively, our data is providing a coherent understanding of how 
subde differences between the KIM-PTPs contribute to MAPK 
selectivity and, in turn, MAPK signal fidelity. 

Experimental Procedures 

Peptide and protein preparation 

The PTPSLjqm peptide (residues 331-348), STEP KIM (residues 
213-229) and STEP KIMKIS (residues 213-256) were synthesized, 
HPLC purified and verified by mass spectrometry (MS) (>98% 
purity; Biosynthesis, Inc.). Expression and purification of ERK2 
and PTPSLkimkis (residues 332-373) was carried out as 
previously described [14,15]. PTPSL (residues 332-655, human), 
STEP (residues 213-539, mouse) and PS chimeril (PTPSL residues 
332-373 with STEP residues 257-539; gene synthesized by DNA 
2.0) were expressed and purified using the protocol previously 
described in detail for STEP CAT (residues 244-539) [19]. The 
DUSP16 MAPK binding domain (residues 5-138, human) was 
expressed and purified as previously described [18]. All purified 
proteins were stored at —80 °C until use at which point they were 
thawed and exchanged into the appropriate buffer using size 
exclusion chromatography (SEC; Superdex 75 26/60; GE 
Healthcare). The final protein buffer was dependent on the 
experiment performed. For ITC experiments, the final SEC 
purification was performed in ITC buffer (10 mM Tris pH 7.5, 
150 mM NaCl, 0.1 mM EDTA, 0.5 mM TCEP); for SAXS 
experiments, the final SEC purification was performed in SAXS 
buffer (50 mM HEPES pH 6.8, 150 mM NaCl, 5 mM 
dithiothreitol). PTPSL K im, STEPiqm and STEPktmkis peptides 
were solubilized in ITC buffer prior to ITC measurements. 

Isothermal Titration Calorimetry 

ITC experiments were performed at 25 °C using a VP-ITC 
microcalorimeter (Microcal Inc.). Titrant (10 uL per injection) was 
injected into the sample cell over a period of 20 seconds with a 250 
second interval between titrations to allow for complete equilibra- 
tion and baseline recovery. 28 injections were delivered during 
each experiment and the solution in the sample cell was stirred at 
307 rpm to ensure rapid mixing. To determine the thermody- 
namic parameters (AH, AS, AG) and binding constants (Kj of the 
ERK2:PTPSL K im, ERK2:PTPSL, ERK2:STEP and 
ERK2:PSChimera interactions, PTPSL K im, PTPSL, STEP and 
PSChimera were titrated into ERK2. To determine the same 
parameters for the ERK2:STEP K im, ERK2:PTPSL K imkis and 
ERK2:STEP K imkis interaction, ERK2 was titrated into the 
peptides. Titrants were chosen based on the amount of purified 
protein/ peptide available for the experiment. For all titrations the 
concentrations were slighdy varied to ensure independence of the 



K d from the protein/peptide concentrations. Data were analyzed 
with a one-site binding model assuming a binding stoichiometry of 
1:1 using Origin 7.0 software. A nonlinear least-squares algorithm 
and the titrant and sample cell concentrations were used to fit the 
heat flow per injection to an equilibrium binding equation, 
providing values of the stoichiometry («), change in enthalpy (AH) 
and binding constant (K). All data were repeated in triplicate. 

Small angle X-ray scattering 

To form the ERK2:STEP resting-state complex, purified, 
unphosphorylated ERK2 was incubated in a 1:1 molar ratio with 
purified STEP 2 13-539 at room temperature for 30 minutes and 
purified using SEC (Superdex 75 26/60 [GE Healthcare]; SAXS 
buffer). Data for the ERK2:STEP complex were recorded on 
samples at 0.7, 1.2 and 1.7 mg/ml. All samples were prepared 
within 48 hours of data acquisition and stored on ice at 4 °C. All 
samples were filtered (0.02 uM filter, Whatman) immediately prior 
to data collection. All data was recorded at Beamline X9A at the 
National Synchrotron Light Source (NSLS) using a Dectris Pilatus 
300k (3.4 m distance from the sample for SAXS) and a Photonic 
Science CCD (0.47 m from the sample for WAXS) detector. 20 ul 
of sample was continuously flowed through a 1 mm diameter 
capillary and exposed to an X-ray beam for 30 seconds. 
Normalization for beam intensity, buffer subtraction and merging 
of the data from both detectors were carried out using PRIMUS 
[20]. A Guinier approximation, I(q) = I(0)exp(-q 2 R g 2 /3), where a 
plot of ln[I(q)] and q is linear for q<1.3/R g (q is the scattering 
vector q = 47tsin(6)/X; X, wavelength of X-rays and 26 is the angle 
between the incident beam and the scattering X-rays) was 
performed on at least four independent scattering trials and 
averaged to determine the radius of gyration (R g ). The linearity of 
the Guinier region and the intensity at zero scattering angle, 1(0), 
were used to validate that all samples were monodisperse in 
solution. I(0)/c, where c is concentration, was consistent for all 
measurements. GNOM [21] was used to determine the pair- 
distribution function, P(r), for the ERK2:STEP resting-state 
complex using data collected on the 1.7 mg/ml sample. Twenty- 
four envelopes were generated using GASBOR [22] and were 
aligned and averaged using the DAMAVER program suite [23]. 

Results 

KIM-PTPs bind ERK2 more tightly than p38a 

We used ITC measurements to determine the binding affinity of 
STEP and PTPSL for unphosphorylated ERK2 (Table 1; 
Figure SI). ITC showed that STEP binds ERK2 ~12-times 
more tightly than p38a. Indeed, the ERK2:STEP complex is the 
tightest binding complex of all MAPKKIM-PTP complexes 
tested. This complex has a ~3-times lower K d (tighter binding) 
than the next tightest complex, ERK2:PTPSL, and a nearly 70- 
fold greater affinity than the weakest complex, p38a:HePTP [15]. 
Like STEP, PTPSL binds ERK2 more tightly (~3. 5-fold) than 
p38ot. We also tested peptides that include the kinase specificity 
sequence (KIS) regions, as we [14,15] and others [12] have shown 
that KIS residues can interact directiy with ERK2 and p38a. Each 
STEP (STEPkim, STEPkimkis, STEP; Figure 1A) and PTPSL 
(PTPSLkim, PTPSLkimkis, PTPSL; Figure 1A) domain has a 
higher binding affinity for ERK2 than for p38a [14]. This 
increased affinity for ERK2 was also observed for each domain of 
HePTP [15,17]. Interestingly, the tightest p38ot:KIM-PTP com- 
plex, p38a:PTPSL (K d ~380 nM), is equivalent to the weakest 
ERK2:KIM-PTP complex, ERK2: HePTP (K d ~370 nM) [15,17]. 
Therefore, when compared to p38a, ERK2 has a higher affinity 
for each KIM-PTP and KIM-PTP domain (KIM, KIMKIS, 
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Figure 1. The ERK2:STEP resting-state complex. (A) Constructs used in this study; (B) SAXS data (l(q) vs q) of the ERK2:STEP resting-state 
complex (black squares); error bars (grey lines). Error bars show the experimental error based on circular averaging of the 2D solution scattering data; 
theoretical scattering curve from calculated ab initio molecular envelope (red); inset, Guinier plots for samples at 1.0 mg/ml and 1.7 mg/ml; (C) The 
ERK2:STEP ab initio molecular envelope in two views rotated by 90 with the dimensions of the envelope. 
doi:1 0.1 371 /journal.pone.0091 934.g001 



KIM-PTP), showing that the increased interaction strength is 
clearly driven by the MAPK rather than the regulatory protein. 

We also generated a chimeric protein in which we fused the 
PTPSL KIM kis to STEP CAT (PS chimera ; Figure 1A) [14]. This was 
done in order to test the contribution of the STEP CAX to the 
MAPKiprotein interaction, as we previously showed that the 
STEPcat plays a role in binding p38a [14]. The PS chime . ra binds 
marginally more tightly to p38a than any of the WT-KJM-PTPs, 
but binds to ERK2 with about equal binding strength as STEP 
(Table 1) [14]. 

The ERK2:STEP resting-state complex is compact 

ITC measurements showed that KIM-PTPs interact more 
strongly with ERK2 than with p38oc. In order to identify if these 
increased interaction strengths are correlated with structural 
changes we used SAXS to determine the overall shape of these 
complexes in solution. We recently showed that the ERK2:HePTP 
resting-state complex is highly extended in solution with a radius of 
gyration (R g ) of 33.3 A and a maximum dimension (D max ) of 
110 A [16]. Here, we used SAXS to determine the size and shape 
of the ERK2:STEP resting-state complex (Figure IB). The 
ERK2:STEP complex was purified by SEC and SAXS data were 
recorded at Beamline X9A at the National Synchrotron Light 
Source (Brookhaven National Laboratories). Guinier analysis of six 
independent SAXS samples was used to calculate an Rg of 



30.6±0.4 A, showing that the R g of the ERK2:STEP resting-state 
complex is —3.0 A smaller than that of the ERK2:HePTP resting- 
state complex (Table 2; Figure 1C; HePTP 324 residues; STEP 
326 residues) [16]. To analyze this difference in more detail, we 
determined the ab initio molecular envelope of the ERK2:STEP 
resting-state complex (Figure 1C). A D max of 95 A was 
determined by analysis of the P(r) function, which is ~ 1 5 A 
shorter than that of the ERK2:HePTP complex (Figure 2) [16]. 
Thus, like the p38a:STEP resting-state complex [14], the ERK2:- 
STEP resting-state complex is compact (Figure 2), and suggests 
that the STEP CA t interacts directly with ERK2. Consistent with 
this hypothesis, STEP binds ERK2 1 1 -times more tightly than the 
STEPkim peptide (Table 1). Based on our ITC data, both the 
STEPkis and the STEP C at contribute to the binding between 
ERK2 and STEP, as the STEPkimkis binds to ERK2 ~4.5-times 
more tightly than the STEP KIM and ~2.5-times less tightly than 
full-length STEP (Table 1). 

The calculated R g and D miLX of the ERK2:STEP resting-state 
complex are highly similar to those determined for the p38a:STEP 
resting-state complex (Figure 2). Using NMR spectroscopy, we 
previously determined that STEP CAT interacts directiy with p38a; 
specifically, helix 0C-2' of STEP binds between the p38ot activation 
loop and the MAPK-insert [14]. Our SAXS and ITC data suggest 
that STEPcat also interacts with ERK2. Although there are a few 
structural differences between ERK2 and p38a in this region (e.g., 
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Table 1. Thermodynamic and dissociation constants for 



ERK2:KIM-PTP and DUSP16 MKB domains derived from ITC 
experiments at 25 °C. 




Complex 


K d InM) 


AH (kcal mol 


TAS (kcal mol 


ERK2: STEPkim 


546±56 


-27.9±0.4 


-19.4 ±0.3 


ERK2: STEP K | MK | S 


123±28 


-29.2±0.8 


-19.7 ±0.9 


ERK2: STEP 


48 ±20 


-29.5±2.5 


-19.7±2.0 


ERK2: PTPSL KIM 


463 ±24 


-23.5±1.1 


-14.9±1.1 


ERK2: PTPSL KIMK , S 


825±190 


-21.6±0.8 


-13.3±0.9 


ERK2: PTPSL 


108±36 


-22.3±2.2 


-12.8±2.1 


ERK2: PS c himera 


47±6 


-18.5 ±0.8 


-8.5 ±0.7 


ERK2: DUSP16 MKBD 


1303±215 


-4.8 ±0.2 


-3.2±0.3 





All experiments were performed in triplicate and values reported are the 
experimental average and standard deviation. 
doi:1 0.1 371 /journal.pone.0091 934.W01 



the activation loop is six residues longer in ERK2), nearly all 
residues that that we showed participate in the p38a:STEP CAT 
interaction are conserved between ERK2 and p38a [14]. Taken 
together, it is very likely STEPcat interacts with ERK2 via a 
mechanism very similar to that observed for p38a. 

The interaction of PTPSL with ERK2 and p38oc is not 
conserved 

While HePTP interacts in an extended fashion with ERK2 and 
p38ot, our data shows that STEP associates with these MAPKs to 
form more compact, globular complexes. Recently, Balasu et al. 
used cross-linking, mass-spectrometry (MS) and molecular mod- 
eling to generate a structural model for the ERK2:PTPSL resting- 
state complex, in which they predict that the complex is compact, 
with ERK2 and PTPSL associating in a side-by-side manner [24] . 
In this model, PTPSL CAT binds ERK2 in an orientation in which 
the catalytic site is positioned near the ERK2 phosphorylation 
loop. This is different from the p38a:PTPSL complex, which we 
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Figure 2. Comparison of the P(r) functions of the MAPK:KIM- 
PTP resting state complexes. p38ccHePTP, maroon; p38oe:PTPSL, red; 
p38oc:STEP, dark blue; ERK2:HePTP : pink; ERK2:STEP, light blue (all 
p38ocKIM-PTP and ERK2:HePTP data were previously published 
[14,15,17], but are shown here to provide a better comparison of the 
P(r) functions). 

doi:1 0.1 371 /journal.pone.0091 934.g002 



Table 2. SAXS analysis of the ERK2:STEP complex. 



ERK2:STEP 



Gulnier approximation 


Rg (A) 


30.6 ±0.5 


P(r) function calculation 


Q-range (A ') 


0.013-0.303 


R g (A) 


31.1 


Dmax (A) 


95 


Structure modeling 


x 2 


2.3±0.2 


NSD 


1.31 ±0.05 



doi:1 0.1 371 /journal.pone.0091 934.t002 

showed is extended in solution [14]. Despite exhaustive attempts 
to collect SAXS data on the ERK2:PTPSL complex, the samples 
persistently exhibited the typical hallmarks of aggregation and 
thus, could not be used for further analysis. However, using ITC, 
we show that full-length PTPSL binds ERK2 ~4-times more 
tightly than the PTPSL KJM peptide (Table 1). This additional 
binding strength is not a result from the interaction of PTPSLkis 
residues with ERK2, as the PTPSL K imkis binds less tighdy than 
the PTPSLkim- Therefore, it is the interactions between 
PTPSLcat an d ERK2 that lead to the increased binding strength 
observed for the ERK2:PTPSL complex. Thus, our ITC data 
support the model Balasu et al. proposed, in which PTPSL CAT 
directly interacts with ERK2 [24] . 

Unlike the KIM-PTPs, DUSP16 binds p38a and ERK2 with 
similar affinities 

The KIM domains of KIM-PTPs are intrinsically disordered in 
solution [14,15]. Thus, they only adopt a single conformation 
when they bind their respective MAPKs. This is different from the 
MAPK-binding domains (MKBD) of MAPK phosphatases (MKP; 
dual specificity phosphatases; DUSPs) [25,26]. In MKPs, the KIM 
sequence is part of a stable, folded domain and thus adopts a 
single, ordered conformation in both the unbound and MAPK- 
bound states [27,28]. Using a combination of biophysical studies, 
we recently showed how the DUSP16 MAPK-binding domain 
interacts with p38ot [18] . Our NMR data showed that the MAPK- 
binding domain of DUSP16 interacts at multiple positions with 
p38ot. While some of them overlap with those used for KIM-PTPs 
(e.g. interaction with <I> B and Y,. hg ) others are distinct (helix 0t4) 
and based on our results important for the interaction of MAPK- 
binding domain of DUSPs with MAPKs. The ITC measurements 
reported a Kj of ~900 nM. In order to determine if DUSP16, like 
the KIM-PTPs, discriminates between p38a and ERK2, we 
carried out similar ITC measurements between the DUSP16 
MKBD and ERK2. ITC reported a K d of -1300 nM (Table 1; 
Figure SI), which is only ~ 1.4-fold weaker than the interaction 
with p38a. Thus, while the KIM-PTPs clearly exhibit an increased 
affinity for ERK2 over p38ot, DUSP16 binds these MAPKs with 
similar affinities. 

Discussion 

Over the last years it has become apparent that the KIM-motifs 
from most regulatory proteins readily discriminate between JNK 
and p38/ERK, however they discriminate much more poorly 
between p38 and ERK. For example, in a recent study, JNK- 
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Figure 3. Regulatory protein-MAPK interaction strengths. The binding strengths of each KIM-PTP (KIM peptides and KIM-containing MKBDs 
(in orange font), light blue; KIMKIS peptides, light purple; proteins, light pink) to p38a and ERK2 represented from weakest to tightest (largest K d to 
smallest KJ binding affinity; '*' indicates values previously published [14,15,17,18]. KIM-PTP:MAPK complexes determined to be extended using SAXS 
written in normal text; KIM-PTP:MAPK complexes determined to be compact written in italics; *, compact nature of the ERK2-PTPSL complex 
determined by Balasu ef al. [24] and corroborated by ITC measurements (Table 1). 
doi:1 0.1 371 /journal.pone.0091 934.g003 



specific KIM peptides were only specific to JNK, with no 
detectable binding of these JNK-specific peptides by p38/ERK. 
In contrast, KIM peptides specific for either p38 or ERK2 
exhibited discrimination factors of less than 4-fold between p38 
and ERK2 [7]. Here, we show that while a subset of p38/ERK2 
regulators, the KIM-PTPs, interact productively with both p38ct 
and ERK2, their affinities for ERK2 are higher (Figure 3). 
Specifically, the KIM peptides of the KIM-PTPs exhibit 
differences in affinity for ERK2 versus p38ot between 4.0 (HePTP) 
and 6.5 (STEP). Furthermore, for the interaction between MAPKs 
and the KIM-PTP proteins (i.e., not just the KIM-peptides), the 
fold differences in affinity become much larger, to as high as 1 2 for 
STEP. This shows that residues outside the KIM contribute to the 
overall interaction strength of MAPKs for their regulatory proteins 
(Figure 3). It also demonstrates that the increased interaction 
strength is clearly driven by the MAPK rather than the regulatory 
protein. The amino acid sequence of p38a is ~45% identical to 
ERK2 and multiple residues that line the KIM hydrophobic 
pocket, especially around 0 B h^S^ER^: G110E/A111T/ 
N115K) and the CD site (p38 GD ERK2: D124N/Q128C/ 
E160T/D161T) are not conserved between the two MAPKs. 
Thus, these sequence differences are the likely reason for the 
reduced affinity of the regulatory proteins for p38a versus ERK2. 
In contrast, and more consistent with the previous study, the 
interaction of the KIM in DUSP 1 6 exhibits almost no difference 
in affinity for p38a and ERK2, with a difference in K<j of only 
~ 1.4-fold. Together, our data shows that while KIM peptides 



typically poorly discriminate between ERK2 and p38(X (as 
compared to JNK), significant discrimination for these two 
MAPKs can be seen for at least a subset of regulatory proteins, 
the KIM-PTPs. Collectively, these studies are providing a more 
coherent understanding of how differences in the family of KIM- 
PTPs with ERK2 or p38a contribute to MAPK selectively and, in 
turn, MAPK signal fidelity. 

Supporting Information 

Figure SI Isothermal titration calorimetry data. Raw 

isothermal titration calorimetry data (upper panels) and derived 
binding isotherm plotted versus the molar ratio of titrant fit using a 
one-site model (lower panels) for ERK2 with (A) PTPSLkim, (B) 
PTPSLkimkis, (C) PTPSL, (D) DUSP 16 MKBD, (E) STEP K im, 
(F) STEPkjmkis, (G) STEP and (H) PS ohimera . 
(TIF) 
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